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Abstract 

The electron and muon number violating muonium-antimuonium oscillation process can 

proceed provided neutrinos have non-zero masses and mix among the various generations. 

Modifying the Standard Model only by the inclusion of singlet right handed neutrino fields 
and allowing for general neutrino masses and mixings, the leading order matrix element 
contributing to this process is computed. For the particularly interesting case where the 
neutrino masses are generated by a see-saw mechanism with a very large Majorana mass 

Mr >> Mw, it is found that both the very light and very heavy Majorana neutrinos each 
give comparable contributions to the oscillation time scale proportional to M R . Present 

experimental limits set by the non-observation of the oscillation process sets a lower limit 

on Mr of roughly of order 10 4 GeV. 
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One of the most striking examples of a purely quantum mechanical phenomenon which 
occurs in a wide range of physical systems involves the time dependent oscillation between 
two distinct levels or particle species. Examples of such systems are quite varied and range 
from the text book example of a particle moving in a double well potential of the ammonia 
molecule to oscillations in the neutral K° — K° and B° — B° meson systems [1]- [3]. The later 
processes arise as a consequence of the fact that the left handed quarks which participate 
in the weak interaction are linear combinations of the mass diagonal quark states. Thus the 
weak interaction currents coupling to the W ± vector bosons when written in terms of the 
quark mass eigenstates contain the left handed charge 2/3 quarks and a mixture of the left 
handed charge -1/3 quarks. Thus the neutral meson can convert into its antiparticle via the 
exchange of two W vector bosons which arises in fourth order in the weak interaction. 

Of the various fermions appearing in the Standard Model, the only ones which are 
electrically neutral are the neutrinos. Left-handed neutrinos are components of SU(2)l 
doublets along with their charged leptonic partners and experience only the weak interaction, 
while any right-handed neutrinos are completely neutral under the Standard Model gauge 
group. During the past several years, there has been mounting experimental evidence [4]- [7] 
that the neutrinos involved in the weak interactions are in fact linear combinations of nonzero 
mass neutrino eigenstates in a somewhat analogous fashion to the quarks. Thus the various 
weak interaction neutrinos exhibit a mixing phenomenon and, being electrically neutral, can 
also oscillate from one species into another. The size and nature of the neutrino mass and 
the associated mixing is still an open question subject to experimental determination and 
theoretical speculation [8]. 

Muonium (M) [9] is the Coulombic bound state of the electron and the antimuon (e~/i + ), 
while antimuonium (M) is the Coulombic bound state of the positron and the muon (e + /i~). 
In order for there to be a nontrivial mixing between these two states, the individual electron 
and muon number conservation must be violated. Such a situation results provided the 
neutrinos are massive particles which mix amongst the various generations. Modifying the 
Standard Model only by the inclusion of singlet right handed neutrinos and allowing for a 
general mass matrix for the neutrinos, this criterion can be met and the e~ ji + and e + fi" 
states can indeed mix [10]. 

The lowest order Feynman diagrams producing muonium and antimuonium mixing are 
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displayed in Figure 1. 
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Figure 1: Feynman graphs contributing to the muonium-antimuonium mixing. 

They contain two neutrinos and two W bosons in the intermediate state. The first two graphs 
(a and b) contributing to the muonium-antimuonium mixing bears a striking resemblance to 
the quark box diagram underlying the mixing in the neutral kaon and B meson systems [2] . 
For the M — M case, however, one is free from the complications of the strong interactions 
since both M and M are simple nonrelativistic Coulombic bound states composed of leptonic 
constituents. Thus the mixing of the M — M system can be unambiguously calculated. On 
the other hand, the later two graphical contributions (c and d) arise since, in general, the 
neutrino mass eigenstates are Majorana (self-conjugate). These graphs have no counterpart 
in the quark case. As it turns out, however, these graphs cancel against each other in 
the calculation of the effective Lagrangian and thus do not contribute to the muonium- 
antimuonium oscillation process. 
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The Standard Model leptonic charged current interaction is given by 

4f = -J—^fw+v + jffw-») , (1) 

SITluw 

where the (purely left-handed) leptonic charged current is 

^-T = (4T) t = E^7^ffi- (2) 

a=l 

Here and v^l are respectively the charged lepton and its associated neutrino partner 
of the SU(2) L doublet. The superscript zero indicates weak interaction eigenstates. A 
generic left handed fermion field is ip L — ^(1 — 75 )ip, while the right handed fermion field 
is ipR — |(1 + 75) ip. The fields participating in the weak interaction, however, need not be 
mass diagonal. The mass term for the charged leptons arises from their Yukawa couplings 
with the scalar doublet. After spontaneous symmetry breaking, the mass term takes the 
form 

£"mass = ~ E l^Ra m ab^Lb + ^La m ba^Rb] > (3) 
a,b=l 

where m e is a 3 x 3 mass matrix. To diagonalize this matrix, one performs a biunitary 
transformation 

m f = A R m^ diag (A L ) t , (4) 

where A R and A L are 3x3 unitary matrices and m diag is a diagonal 3x3 matrix whose 
entries are the charged lepton masses, rri£ a . To implement this basis change, the charged 
lepton fields participating in the weak interaction are rewritten in terms of the mass diagonal 
fields (denoted without the superscript) as 

'2 = E A L ab £ Lb ; 42 = E A* in . (5) 

a=l a=l 

So doing the mass term reads 
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C-Lass = - E m la\tRa£ha + ha^Ra] • (6) 
a=l 

A general neutrino mass term takes the form 
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L v ""a L l^La L " m ab U Lb + ^La Um ba ^Lb J _ Z^rffi m ia V La + V La m ia V Ri \ 
Z a,6=l a=l i=4 
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2 E[4?^g^ T + 4? T ^^]. (7) 



i,J=4 



Here m L ,m D ,m R are Lorentz and SU(2) L x {7(1) singlet 3x3 complex matrices. m D 
preserves lepton number, but violates SU(2) L x {7(1) since it connects an SU(2) L doublet 
with a singlet and thus can only be generated after spontaneous symmetry breaking. This 
lepton number conserving mass is often referred to as a Dirac mass. Both m L and m R violate 
the individual lepton numbers by two units and are referred to as Majorana masses. The 
presence of m L also requires SU(2)l x {7(1) symmetry breaking since it connects left handed 
neutrinos which are parts of an SU(2)l doublet. In the Standard Model its generation 
requires a mass dimension five operator. On the other hand, since m R connects SU{2) L x 
{7(1) singlet right handed neutrino fields, it can appear even in the absence of SU(2) L x {7(1) 
symmetry breaking. 

The matrix C appearing in the Majorana neutrino mass terms is there to preserve Lorentz 
invariance. This matrix satisfies 

CC* = 1 ; C = C T ■ C^C' 1 = - 7 „ . (8) 

The specific form for C is representation dependent. In the Dirac representation, C = 2727°. 
The charge conjugate field i/j c is defined as 

^ c ( x ) = C^ T {x) ; if = i) T C . (9) 

Using the charge conjugated fields, one can write 

1, 



•M = §[*M + ("Sw^ri do) 

and thus it follows that the neutrino mass term can be cast in the compact form as 



For 3 generations of neutrinos, the six mass eigenvalues, m u A, are obtained from the diago- 
nalization of the 6x6 matrix 

/ m L (m D ) T \ 

= D R ■ (12) 
\ m m / 

Since 1VP is symmetric, it can be diagonalized using a unitary transformation 

M^ iag = (U)tM"U , (13) 
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where U is a 6 x 6 unitary matrix. This diagonalization is implemented via the basis change 
on the original neutrino fields organized as the 6 dimensional column vector 

to the new neutrino fields (without superscripts) defined as 

NW=UN L ; N { V = UN R , (15) 

where 

The neutrino mass term takes the form 

I 6 6 

^mass = ^E m vA \v T A Cv A + va.Cv t a \ = - E m^PJ^ , (17) 

Z A=l A=l 

where are the Majorana neutrino masses. Note that these mass diagonal fields are 
Majorana (self-conjugate) fields satisfying v A = (N L + N R ) A = v c A ; A = 1, 6. 

As a consequence of the diagonalization of the neutrino mass matrix, the charged current 
when written in terms of the mass diagonal fields takes the form 

J-Z = (4T) f = EE h~a^V aA V LA , (18) 
a=l A=l 

where 

V aA = jz{A L l ) ab U bA . (19) 

6=1 

Thus the neutrino fields appearing in the charged weak current are given by 
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v 



(0) 
La 



E UaAVLA , (20) 
A=l 

where v LA is the A th component of N L . Note that since both A L and U are unitary matrices, 
it follows that 

E v aA v; A = 5 ab . (21) 

A=l 

Since a nonzero Majorana mass matrix m R does not require SU{2)l x U{\) symmetry 
breaking, it is naturally characterized by a much larger scale, M R , than the matrices m D 



and m L whose nontrivial values do require SU(2)l x U(l) symmetry breaking. In fact, since 
the presence of m L requires a mass dimension five operator, it is naturally much smaller 
than both Mr and the elements of m D which are expected to be somewhere of the order of 
the charged lepton mass to the W mass. Thus for simplicity, one can set m L to zero, while 
taking the elements of m D , characterized by a scale mo, to be much less than Mr, the scale 
of the elements of m R . So doing, one finds on diagonalization of the 6x6 neutrino mass 
matrix that 3 of the eigenvalues are crudely given by 

m va ^-^«m D ■ a = 1,2,3, (22) 

Mr 

while the other 3 eigenvalues are roughly 

m ui ~ M R ; i = 4,5,6 . (23) 

This constitutes the so called see-saw mechanism [11] and provides a natural explanation of 
the smallness of the 3 light neutrino masses. Moreover, the elements of the mixing matrix 
are characterized by the Mr mass dependence 

U ab ~0(l) ; o,6=l,2,3 
t%~0(l) ; i,j = 4,5,6 

U*~UZ~0(^) ; a = 1,2,3 ; i = 4,5,6. (24) 

Mr 

Since the charged lepton mixing matrix is independent of M R , one finds that elements of 
the mixing matrix appearing in the charged current has the Mr mass dependence 

V ab ~0(l) ; a,b= 1,2,3 

V ai ~0(^) ; a =1,2,3; i = 4,5,6. (25) 

Mr 

Note that the states v a ; a = 1,2,3, are predominately composed, up to corrections 
of order ^ << 1, of the neutrino fields vfl which are the ones participating in the weak 
interaction. Thus for most practical applications, one can simply neglect the heavy neutrino 
fields which are primarily composed of the fields v^- ; i — 4, 5, 6, which are electroweak 
singlets. As we shall see, however, when considering the case of muonium-antimuonium 
oscillations, the heavy neutrino states will yield contributions of the same order as the light 
neutrino states and their presence cannot be neglected. 
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The evaluation of the graphs of Fig. 1 contributing to the muonium-antimuonium mixing 
is straightforward. Taking the external leg momenta to vanish, the result of the calculation 
can be encapsulated in an effective Lagrangian [12] of the form 

£eff = ^JfWil ~ 75)e][7Z7,(l - 75)e] . (26) 
An explicit unitary gauge calculation yields 

6 

(27) 



r< r<2 ]\/f2 



£ {V, A V} A fS{x A ) + J2 {V»aV} a ){V, b V} b )T{ x A ,x B ) 

A=l A,B=1;A^B 



^2 16tt 2 

where G F ~ 1.16 x 1(T 5 GeV~ 2 is the Fermi scale and x A = ^ ; A = 1, 6. Here 

„. . x 3 -llx 2 + 4x 3x 3 n . . 

= ^tr^r — w^ tn(x) (28) 

is the Inami-Lin [13] function and 

T(x A , x B ) = x A x B ( J(yXj ^ J ( Xb ) j = t(x b , x A ) (29) 

\ x A - X B ) 

with 

J (x) = ^- 8x+ J>en( X ) - ? * . (30) 
K ' 4(1 -x) 2 v ; 4(1 -x) K ' 

In obtaining this result judicious use was made of the unitarity of the charged current mixing 
matrix elements (c.f. Eq.(21)). As a consequence of the fermi statistics, the graphs of figures 
(a) and (b) gave identical contributions to the effective Lagrangian, while the graphs of 
figures (c) and (d) gave canceling contributions. Note that in order for the neutrinos in 
the intermediate state to give a nonvanishing contribution to the effective Lagrangian, they 
must be massive as well as exhibit a nontrivial mixing with both the electron and the muon 
in the charged current. 

Muonium (antimuonium) is a nonrelativistic Coulombic bound state of an electron and 
an anti-muon (positron and muon). The nontrivial mixing between the muonium (\M > ) 
and antimuonium (|M >) states is encapsulated in the effective Lagrangian of Eq. (26) and 
leads to the mass diagonal states given by the linear combinations 

|M±>=-==l^=[(l + e)|M >±(l-e)|M>] (31) 
y^l + H 2 ) 

where 
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with 

< M\-Jd*r£ c{i \M> <M\-Jd*rC c{i \My 
JVI MM ~ I ; — , — — i JVI MM — I ; — , — — • \ 66 > 

y< M\M >< M\M > \J < M\M >< M\M > 

Since the neutrino sector is expected to be CP violating, these will be independent, complex 
matrix elements. If the neutrino sector conserves CP, with \M > and \M > CP conjugate 
states, then M. M jj = -Mj^m an d e = 0. In general, the magnitude of the mass splitting 
between the two mass eigenstates is 



|AM| = 2\R^M Wm M MTi \ . (34) 

Since muonium and antimuonium are linear combinations of the mass diagonal states, an 
initially prepared muonium or antimuonium state will undergo oscillations into one another 
as a function of time. The muonium-antimuonium oscillation time scale, tjj m , is given by 

—!— = \AM\ . (35) 

T MM 

A nonrelativistic reduction of the effective Lagrangian of Eq. (26) produces the local, 
complex effective potential 

V cS (f) = 8^5 3 (r) . (36) 

Taking the muonium (antimuonium) to be in their respective Coulombic ground states, 
0ioo if) = J-^ er r l a ~iM , where ajj M = — is the muonium Bohr radius with m re d = 
~ m e the reduced mass of muonium, it follows that 



m e +m M 

|Re GxfiiJ , , ,„x,o 16]ReG 



— ~ 2 / d 3 r0* loo (f)|Re V cS (r) |0 loo (f) = 16 |Re G " m1 |0 1OO (O)| 2 = - 



MM I 



'MM " n V<Z °^MM 

(37) 



Thus we secure an oscillation time scale 

_L ~ l 3 ReG MM\ m , a , 

T MM n V2 

with ^M- given by Eq. (27). 



(38) 
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The present experimental limit [14] on the non-observation of muonium-antimuonium 
oscillations translates into the bound |Re Gjj m \ < 3.0 x 10~ 5 G,f. This limit can then be 
used to construct a crude lower bound on Mr. For the case when the neutrino masses arise 
from a seesaw mechanism and taking mo to be of order M w , the Mr dependence of Gjj m 
is obtained from Eq.(27) as: 

G 2 M 4 

I^GmmI^-T^ ; A = a = 1,2,3; 5 = 6 = 1,2,3 

R 

Q2 j\ j4 

|ReG MM |~^-^ ; A = i = 4,5,6 ; B = j = 4,5,6 



iReG^I-^^M^) 5 A = a = 1,2,3; 5 = ^ = 4,5,6. (39) 
That is, the contributions to |Re Gjj M \ are of comparable magnitude when the neutrinos 

wrv Her tit m ^, f)( 



propagating in the loop are both either very light m a ~ Oi^f-) or very heavy ~ Mr, 



but are suppressed by an additional factor of ^^£n(j^) when one of the neutrinos is light 
and the other neutrino is heavy. For the two light neutrino intermediate state, there is one 
factor of coming from each of the two light masses, while the mixing matrix elements 
are of order unity. On the other hand, for the case when there are two heavy neutrinos in 
the intermediate state, there is a factor of arising from each of the four mixing matrix 
elements and two factors of the heavy mass, Mr, from the two masses once again resulting 
in a net factor of Thus, even though the heavy Majorana neutrinos constitute but a 
very small amount (C(f^)) of the neutrinos which participate in the weak interactions, 
they still contribute at the same level (in terms of M R dependence) as the light neutrinos 
which are the principle components of the weak interaction neutrino fields. 

Unfortunately, the current experimental bound on |Re Gj^ M \ is not particularly restric- 
tive yielding the modest bound 

Mr > C(10 4 GeV) . (40) 

We stress that this bound is at best an order of magnitude estimate since we are retaining 
only the mass dependence on Mr and neglecting all numerical dependence on the mixing 
angles and CP violating phases in V a A- 

This work was supported in part by the U.S. Department of Energy under grant DE-FG02- 
91ER40681 (Task B). 
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